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, Abstract. Adaptive SPH and N-body simulations were carried out to study the 

' evolution of the equilibrium structure of dark matter halos that result from the grav- 

itational instability and fragmentation of cosmological pancakes. Such halos resemble 
those formed by hierarchical clustering from realistic initial conditions in a CDM uni- 
verse and, therefore, serve as a test-bed model for studying halo dynamics. The dark 
matter density profile is close to the universal halo profile identified previously from 
^) , N-body simulations of structure formation in CDM, with a total mass and concentra- 

' tion parameter which grow linearly with scale factor a. When gas is included, this 

, concentration parameter is slightly larger than the pure N-body result. We also find 

that the dark matter velocity distribution is less isotropic and more radial than found 
by N-body simulations of CDM. 
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' CDM Simulations vs. Observed Halos N-body simulations of structure for- 

■ mation from Gaussian-random-noise density fluctuations in a cold dark matter 
(CDM) universe have revealed that dark matter halos possess a universal density 
profile that diverges as r~'^ near the center, with 1 < 7 < 2 [6] [7]. There exists 

^ ' a discrepancy between these singular density profiles found in N-body simulations 

■ and current observations of the rotation curves of nearby dwarf galaxies [6] and 
of strong lensing of background galaxies by the galaxy cluster CL0024-I-1654 [10] 
[8], which suggest that dark- matter dominated halos of all scales have fiat density 
cores, instead. 

Halo Formation by Pancake Instability and Fragmentation The model we 
use to examine the formation of dark-matter dominated halos is that of cosmological 
pancake instability and fragmentation, previously discussed in detail by [12]. Halos 
formed by such a pancake instability have density profiles very similar to those 
formed hierarchically in CDM models (e.g. NEW profile [7]), providing a convenient 
alternative to more complicated simulations with more realistic initial conditions 
[1] [5] [11]. The ASPH/P'^M simulations considered here were described by [1]; 
this paper extends that analysis to evolutionary trends in the dark matter halo 
structure. 
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FIGURE 1. (left) Dark matter density field at a/ac ~ 3. (middle) Density profile of the dark 
matter halo as simulated without gas at four different scale factors, a/oc = 3 (solid), 4 (dotted), 5 
(short dash), and 7 (long dash). Shown above are fractional deviations [pnfw — p)/ Pnfw from 
best-fit NFW profile for each epoch, (right) Same as middle, but for DM halo simulated with 
gas+DM. 

Main Results: 

• For a/oc between 3 and 7, the halo can be fit by an NFW profile, with mass 
within r2oo growing linearly with scale factor a, when simulated either with 
or without gas: M2qq{x) ~ 0.07a;, where x = a/ttc, and Oc is the scale factor 
at primary pancake collapse (see Figs. 1 & 2). This mass evolution resembles 
that of self-similar spherical infall [2], despite the anisotropy associated with 
pancake collapse and filamentation and periodic boundary conditions. 

• After a/a^. = 3, the concentration parameter c^fw = '"2oo/'"s) determined by 
best-fitting an NFW density profile [7] to our simulation halos, grows roughly 
linearly with scale factor a: cnfw{x) ~ 1.33a; — 0.18 (without gas), cnfw{x) ~ 
1.49a; — 0.37 (with gas) (i.e. the linear slope is steeper in the case with gas 
included). Fluctuations in cnfw around this trend are smaller when gas is 
included. This evolution we find for c^pw is reminiscent of that reported for 
halos in CDM N-body simulations [3]. However, the latter applies to halos of 
a given mass which are observed at different epochs, and, therefore, refiects 
the statistical correlation of halo mass with collapse epoch in the CDM model, 
while our result follows an individual halo. 

• The anisotropy parameter (3 = 1 — {Vf)/{2{v^)), where Vt{vr) are tangential 
(radial) velocities, is shown in Figure 2. Pancake halos are somewhat radially 
biased, with (3 > 0.6, about twice the value reported for halos in CDM N-body 
simulations [4] [9]. With no gas included, the average anisotropy in the halo 
does not change very much with time, while the inclusion of gas leads to a 
slight drop after a/uc — 5. 




FIGURE 2. (left) Evolution of halo dark matter integrated mass Mx as simulated with (dotted) 
and without (solid) gas, within spheres of average ovcrdcnsity X = {p)/p (in computational units, 
where Mi,ox = ^pP = !)• (middle) Evolution of halo concentration parameter for the dark matter 
halo as simulated with (dotted) and without (solid) gas. (right) Anisotropy parameter /3 averaged 
over all dark matter halo particles within a sphere of average overdensity 200, as simulated with 
(dotted) and without (solid) gas. 

This work was supported by NASA ATP grants NAG5-7363 and NAG5-7821, 
NSF grant ASC-9504046, and Texas Advanced Research Program grant 3658-0624- 
1999. 



REFERENCES 

1. Alvarez, M., Shapiro, P.R., & Martel, H., RevMexAA (S.C.), in press (2001) (astro- 
ph/0006203) 

2. Bertschinger, E., ApJS, 58, 39 (1985) 

3. Bullock, J.S. et ai, MNRAS, 321, 559 (2001) 

4. Eke, V.R., Navarro, J.F., Frank, C.S., ApJ, 503, 569 (1998) 

5. Martel, H., Shapiro, P.R., &: Valinia, A., in preparation (2001) 

6. Moore, B., Quinn, T., Governato, F., Stadel, J. & Lake, G., MNRAS, 310, 1147 
(1999) 

7. Navarro, J.F., Frenk, C.S., and White, S.D.M., ApJ, 490, 493 (1997) 

8. Shapiro, P.R., Iliev, I.T., ApJL, 542, 1 (2000) 

9. Thomas, P.A. et al, MNRAS, 296, 1061 (1998) 

10. Tyson, J. A., Kochanski, G.P., and deh'Antonio, LP., ApJL, 498, 107 (1998) 

11. Vahnia, A., Ph.D. Thesis, University of Texas (1996) 

12. Vahnia, A., Shapiro, P.R., Martel, H., & Vishniac, E.T., ApJ, 479, 46 (1997) 



